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We present a measurement of the ^ oscillation frequency using 
a flavor tagging method based on correlations of B meson flavor with the 
charge of other particles produced in pp collisions at ^/s = 1.8 TeV. Such 
correlations are expected to arise from h quark hadronization and from B** 
decays. We partially reconstruct B mesons using the semileptonic decays 
B^ l'^D^*^~ X and B^ i^D'^X. From the oscillation frequency, we 
obtain the mass difference between the two B^ mass eigenstates, Am^ = 
0.47llQogg(stat) ± 0.034(syst) h ps~^, and measure the efficiency and purity 
of this flavor tagging method for both charged and neutral B mesons. 

PACS numbers: 13.20.He, 14.40.Nd 

The meson is a bound state of a 6 quark and a d quark. Second order electroweak 
processes in which b ^ d while d b result in the transition of mesons into B^ mesons 
and thus B^ ^ B^ oscillations. The b —>■ d coupling occurs via a virtual top quark, thus the 
frequency of these oscillations is sensitive to the magnitude of the element Vtd of the Cabbibo- 
Kobayashi-Maskawa matrix [|I|. For an initially pure sample of B^ mesons (at t = 0), the 
numbers of B^ and B^ mesons at proper time t are given by 

iV±W = A^+(0)^(l±cosAmdt), (1) 

where r is the lifetime of the 5° meson and Am^ is the mass difference between the mass 
eigenstates of the B^-B^ system. In this paper we determine Am^ by measuring the fre- 



quency of flavor oscillations using partially reconstructed semileptonic decays of 

the B meson to iD^*^X. The data used in this analysis were collected with the Collider 
Detector at Fermilab (CDF), at the Tevatron pp Collider at ^/s = 1.8 TeV, and correspond 
to an integrated luminosity of ~ 110 pb^^. 

To extract Arrid the proper time of the B decay is required, as well as the flavor of the B 
meson at the times of its decay and production. While the B flavor at decay is determined 
by its decay products, the determination of the initial B flavor is experimentally challenging. 
Techniques to determine the initial B flavor in several previous measurements of Arrid ^ 
relied on identifying the flavor from the other b hadron in the event {e.g., using the lepton 
charge from the semileptonic decay of this other hadron), and are thus referred to as Opposite 
Side Tagging (OST). 

It has been suggested that the electric charge of particles produced near a B meson 
could also be used to determine its initial flavor, providing a basis for Same Side Tagging 
(SST). For example, if a 6 quark combines with a u quark to form a B^ meson the remaining 
u quark may combine with a d quark to form a vr" (Throughout this paper, reference to 
a specific particle state implies the charge conjugate state as well). Similarly, if a 6 quark 
hadronizes to form a B^ meson, the associated pion would be a vr^. Another source of 
correlated pions are decays of the orbitally excited {L = 1) B mesons (B**) |Q, i.e., B**^ — >■ 
^(*)+7]-^ or B**^ B^*^^7r^. In a hadron collider experiment with central rapidity coverage 
such as CDF, SST methods are attractive since they are expected to have significantly higher 
efficiency than the OST methods. These methods can also be used for measurements of CP 
violation [||. 

In this paper we extract Arrid by applying an SST method to a sample of events containing 
a lepton and a reconstructed D meson from B decay. To determine the initial flavor of 
the B meson, we select one charged track that we will generically refer to as a "pion", 
and use its charge as a tag. We do not attempt to distinguish the hadronization pions 
from those originating from B** decays. The lepton charge tags the 5-flavor at decay 
time. We classify the B-tt combinations as right-sign (RS: B^ti^ and B^n^) or wrong- 
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sign (WS: B+vr"*" and B^tt'). We form the asymmetry in the RS and WS combinations, 
A{ct) = {Njis{ct) — NiYsict))/ {Njis{ct) + N\Ys{ct)), as a function of the proper decay length 
ct. For mesons, we expect an asymmetry independent of ct: A'^lct) = constant = P+. 
The quantity V+ is called the dilution and it is a direct measure of the SST purity, i.e., 
(1 + I^)/2 is the fraction of correctly tagged events. Due to B^B^ mixing, A{ct) for the 
neutral B mesons will vary as a function of the proper decay length ct. From eq. ([l|) and 
the definition of the B^ asymmetry, it follows that the latter is expected to oscillate as 
A^{ct) = Vq ■ cos(Amdt). Mistags, i.e., incorrect flavor determinations, result in a decrease 
of the oscillation amplitude by the B^ dilution factor "Dq- We measure the asymmetry as a 
function of the proper decay length ct for both B^ and B^ mesons, and fit them with their 
expected time dependence, obtaining Am^, Vq and Ref. discusses several effects that 
can result in Vq ^ and in this analysis we impose no constraints on their relative values. 

The CDF detector is discussed in detail elsewhere , and only the features most relevant 
to this analysis are described here. In CDF, the positive z axis is pointed in the proton 
direction, 6 is the polar angle and is the azimuthal angle. The pseudorapidity, rj, is 
defined as — In [tan(0/2)]. We use two devices inside a 1.4 T solenoid to measure tracks from 
charged particles: the central tracking chamber (CTC) and the silicon vertex detector (SVX). 
The combined CTC+SVX tracking system covers the pseudorapidity interval \ri\ < 1.1, 
and gives a resolution on the transverse momentum with respect to the beam axis, p^, of 
S{Pt)/Pt = ((0.0066)^ + (O.OOOQpx)^)^^^ and a resolution on the track impact parameter, 
do, defined as the distance of closest approach to the beamline, of about (13 + 40/pt) /^ni, 
where is in GeV/c. Electromagnetic (CEM) and hadronic (CHA) calorimeters are located 
outside the solenoid and are surrounded by the central muon chambers (CMU) followed by 
the central upgrade muon chambers (CMP). 

The data were recorded using an inclusive lepton (e and fi) trigger. The threshold 
for the single electron trigger was 8 GeV, where Et = E sin 6, and E is the energy measured 
in the CEM. The single muon trigger required a charged track with > 7.5 GeV/c in the 
CTC with matched track segments in both the CMU and CMP systems. Details of the 
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identification of electrons and muons are described in references and 0. 

We use tlie decay chains B° i^i^D^*''^, witli K^n^ir^, or D*^ -D^tt^ followed 

by decaying to K^n^, A'+tt^tt+tt^, or K^tt^tt^, where tt^ denotes the low-momentum 
(soft) pion from the D*^ decay. We use ui~^D^, with K^tt" , where the is 

required not to form a D*~ candidate with another vr candidate in the event. We reconstruct 
the D meson candidates using tracks in a circle of unit radius in 7^-0 space around the 
lepton. For details, see [|l^. We use only tracks that include SVX information. The decay 
products from the D mesons are required to be significantly displaced from the interaction 
point (primary vertex) of the event. The mass distributions of the four decay signatures 
with fully reconstructed D mesons are shown in Fig. |l]a, b and c, while the distribution of 
Am = m{KT{'K^) — m{K7i) for the signature with D*^ -D°vr~, followed by Z)° K^-k^-k^ 
(the 7r° is not reconstructed) is displayed in Fig. |l]d. 

To select the SST pion, we consider all tracks that are within the circle in rj-ip space 
of radius 0.7 centered around the direction of the B meson, approximated by •p{€) +p{D). 
SST candidate tracks should originate from the B production point (the primary vertex of 
the event), and are therefore required to satisfy do/cTdg < 3, where cTrfg is the uncertainty 
on do- String fragmentation models |TT| indicate that particles produced in the 6-quark 
hadronization chain have low momenta transverse to the direction of the 6-quark momentum. 
We thus select as the tag the track that has the minimum component of momentum, 
orthogonal to the momentum sum of the track and the B meson. We define the tagging 
efficiency, e, as the fraction of B candidates with at least one track satisfying the above 
requirements. We measure e ~ 70%, independent of the decay signature used. On average, 
there are about 2.2 SST candidate tracks per B candidate. 

For each of the five decay signatures, we subdivide the candidates into six bins in proper 
decay length, ct. The measurement of ct begins with finding the D decay vertex. The D 
candidate trajectory is intersected with the lepton track (and the momentum of vr* from the 
D*, if present) to form the B decay vertex. We define L^^ to be, in the plane transverse to 
the beam axis, the distance between the primary and the B vertex projected along the B- 



momentum vector. We estimate from the pr of the visible decay products, = |pt(^) + 
p^{D^*^)\, and the mean of the distribution of the momentum ratio /C = {p^ /p^) obtained 
from a Monte Carlo simulation |12[. The p^ /pl^ distribution has a mean of ~ 85% and RMS 
of ~ 12%. The proper decay length of the B meson is then given by ct = L^Juib lp^)KL [p!3|. 



We measure the asymmetry A^'^\ct)^ in each ct bin by simultaneously fitting the mass 
distributions for candidates tagged with a RS or WS pion. The measured asymmetries are 
shown in Fig. ^. If the i.^D'^ and £~^D^*^^ signatures were pure signals of and B^ decays, 
we could simply extract Arrid using the time-dependence of A^{ct). However, the signatures 
are mixtures of B~^ and B^ decays, and thus A^'^\ct) is a linear combination of the true 
asymmetries A^{ct) and A~^{ct). To extract Arrid, "Po and 1^+, it is necessary to determine the 
sample composition of each i^D^*^ signature, by which we mean the fractions of the i~^D^*^ 
candidates originating from the decays of the B^ and B^ mesons. Because a B~^ is associated 
with a vr", whereas an unmixed B^ is associated with a vr"*", the observed asymmetries are 
reduced by cross-contamination. To obtain the true asymmetries A^ and A'^, we introduce 
sample composition parameters (/**, Py, (.norm, Rf, e(7r*), and tb+Ztbo, defined below) to 
describe this cross-contamination and vary these parameters in a simultaneous fit to all of 
the observed asymmetries. 

Cross-contamination can arise if the soft pion vr" from the D*~ decay is not identified - 
the decay sequence 5° ui^D*^ will be reconstructed as i~^D^, that is, as a B~^ candidate. 
We quantify this effect by the efficiency for reconstructing the soft pion, e{7T^). We determine 
e(7r*) by comparing the fraction, R*, of all i^D^X candidates for which we find a vr^ candi- 
date, with the prediction from the other sample composition parameters. (The final result 
is e(7r=K) = 0.85 ± 0.07). Since Z)*° does not decay to D^n~^, there is no cross-contamination 
from B^ ui^D*^ into the B^ sample. 

Another source of cross- contamination arises from semileptonic B decays involving P- 
wave D** resonances as well as non- resonant D^*^n pairs, which cannot be easily recognized 
and removed from the sample. For example, the decay sequence 5° ui^D**^ , followed 
by D**^ — s> D^TT^^ (where by vr** we denote the pion originating from D** decay) will be 
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reconstructed as £"'"-D°, because of the missed vr"; again, a 5° decay is misclassified as 
a candidate. We quantify this effect by the parameter /**, which is the ratio of the 
branching fraction B{B —>■ uiD**) to the inclusive semileptonic B branching fraction, B{B 
iuX). We use /** = 0.36 ±0.12 0. The fractions of 5+ and B° decays in each decay 
signature are also affected by the relative abundance of the four possible spin-parity D** 
states, some of which decay only to D*7r and others to Dn. We define Py = B{B D** 
D*Tx)/ [B{B — s> D** — s> D*7r) + B{B D** Dti)], which we leave as a free parameter in 
our fit. 

The sample composition also depends on / and /*, the ratios of the branching fractions 
B — > uiD and B uiD* to the inclusive semileptonic B branching fraction. We define 



Rf = f*/f, use Rf = 2.5 ±0.6 HTSf, and assume f + f* + f** = 1. The lifetime ratio tb+/tbq 
is another sample composition parameter: Tb+ 7^ t^o implies B{B^ i^uX) 7^ B{B^ 
i'^h'X), as well as a ct-dependent sample composition. We use tb+/tbo = 1.02 ± 0.05 ||15| . 

The tagging is further complicated when a nf^ from D** decay is present. The nf^ may be 
incorrectly selected as the SST pion, always resulting in a RS correlation. The requirement 
do/<^do < 3 reduces this effect: the vr^,* originates from the B meson (iecay point, whereas the 
appropriate tagging track originates from the B meson production point. The contamination 
of the tagging pions from vr^^ decreases with increasing ct. It is quantified by the parameter 



^, defined as the probability of selecting the vr^^, as the tag in a tagged event in which a tt^ 



was produced. We use a Monte Carlo simulation [|T6[ of b quark fragmentation and decay 
to determine ^ as a function of proper decay time. We then use our data to determine a 
normalization factor for ^, C,norm, from the fraction, R**, of all tagged events in which a Tr*,^ 
was selected as the tag. We measure R** from the distribution of impact parameter, (i^, of 
the SST pion with respect to the B decay point. In this study, the do/ao < 3 requirement 
was removed to increase the number of events. We fit the RS distribution of dB/crdB with 
the sum of the scaled WS distribution and a Gaussian of unit width centered at zero. The 
area of this Gaussian yields the total number of SST pions originating from D** decays. 
We extract C,norm by comparing R** with the prediction from the other sample composition 
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parameters. We measure separate values of R** for iD and iD* signatures, and their ratio 
constrains Py- {(.norm = 0.8 ± 0.2 and Py = 0.3 ± 0.3 are the final results). 

The mass difference Arrid and the dilutions T>q and P+ are determined from a fit to the 
measured asymmetries A^"^\ In the fit, /**, Rj and tb+/tbo vary within their uncertainties, 
while e(7r^,), Py, and C,norm vary within constraints imposed through the measured values of 
R* and R**. The prediction for an observed asymmetry is 

Ak = alA^ + at{~A^)+al*{+l) (2) 



where all information about the sample composition is contained in the coefficients a^' [|13 . 
The second term in describes the cross-contamination occurring with the minus sign 
due to the opposite charge correlation for B^, while the third term describes the effect of 
selecting the vr^ as a tag. 

The result of the fit, overlaid onto the measured asymmetries, is displayed in Fig. ^ 
The oscillation in the neutral B signatures is clearly present, giving Arriii = 0.471 h ps^^. 
The fit has a of 26.5 for 30 degrees of freedom. We assign the uncertainties from the 
sample composition parameters, which are dominated by the uncertainty on /**, to the 
systematic uncertainties on Am^, Vq or 1^+. Other sources of systematic uncertainty are 
much smaller, and are due to uncertainties in the Monte Carlo simulation as well as the 
presence of physics background processes that can mimic the i^D^*^X signatures. The 
systematic uncertainties are given in Table |I[ The final result for the mixing frequency 
is Anifi = 0.47llQ;og8(stat) ± 0.034(syst) h ps^^. We also obtain the following values for 
the neutral and charged meson tagging dilutions: Vq = 0.18 ± 0.03(stat) ± 0.02(syst) and 
V+ = 0.27 ± 0.03(stat) ± 0.02(syst). The fit indicates that ~ 82% of the i+D^X signature 
comes from B~^ decays, while ~ 80% of the i^D'X and ~ 95% of the i^D*~'X originate 
from B^. The B^ component of the i~^D^X signature can be seen as a small anti-oscillation 
in Fig. 0, top. 

In conclusion, we have applied a Same Side Tagging technique to samples of 5° — >■ 
£+/)(*)-X and B+ i+D^X decays in pp collisions. The measurement of the asymmetry 
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between tag-charge and S-flavor as a function of proper time for neutral B mesons results 
in the observation of a time- dependent oscillation ^ B^, with the oscillation frequency 
Amrf = 0.47llo;ogg(stat) ± 0.034(syst) h ps~^, in good agreement with the world average 
Arud = 0.474 ±0.031/;, ps~^ |T^. This establishes the effectiveness of Same Side Tagging for 
the first time in hadronic collisions. 
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The version 5.7 was used. We simulate i?-meson decays using the QQ package [12 



Source 


a(P+) 




a{Amd){h ps ^) 


Composition 


+0.0216 
-0.0131 


+0.0225 
-0.0131 


+0.0295 
-0.0310 


Other sources 


±0.0068 


±0.0074 


±0.0146 


Total 


+0.0226 


+0.0237 


+0.0329 


-0.0147 


-0.0150 


-0.0343 



TABLE I. The breakdown of the systematic uncertainties. 
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FIG. 1. Invariant mass distributions of the five B t^D^*^X signatures used in the analysis. 
The sohd histograms correspond to the "right" lepton-kaon charge correlation and the 

dashed histograms the "wrong" one {t'^K~). The solid lines are the fits to the £+Er+ distributions, 
and the resulting number of signal events is given on each figure, (a) K+7r~ mass in ^'^D^X 
K+TT-). (b) K+TT-TT- mass in t+D-X (D" K+tt-tt-). (c) candidate mass in 
£+D*-X {D*- D^TT-, L>o ^ K+TT- and L»° -> K+n-n+n-). (d) K+tt-tt- - i^+vr" mass 
difference in £+D^*^-X {D*- D^tt' , D° K+tt-tt^). The e+K+ background shape (dotted 
line) was determined from the fit (lower solid curve) of the l+K" distribution. 
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FIG. 2. Measured asymmetries as a function of the proper decay length, ct, for the decay signa- 
tures: £^D^ (dominated by B^), £^D~ and the sum of all three £^D*~ (dominated by B^). We 
fit the three £'^D*~ signature separately, but combine them for display purposes. The dashed line 
is the result of the fit. 
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